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ABSTRACT 



Aims. We present a high-spatial and high-spectral resolution observation of the well-studied K giant Aldebaran with AMBER at the 
Very Large Telescope Interferometer (VLTI). Our aim is to spatially resolve the outer atmosphere (so-called MOLsphere) in individual 
CO first overtone lines and derive its physical properties, which are important for understanding the mass-loss mechanism in normal 
(i.e., non-Mira) K-M giants. 

Methods. Aldebaran was observed between 2.28 and 2.31 (im with a projected baseline length of 10.4 m and a spectral resolution of 
12000. 

Results. The uniform-disk diameter observed in the CO first overtone lines is 20-35% larger than is measured in the continuum. We 
have also detected a signature of inhomogeneities in the CO-line-forming region on a spatial scale of ~45 mas, which is more than 
twice as large as the angular diameter of the star itself. While the MARCS photospheric model reproduces the observed spectrum 
well, the angular size in the CO lines predicted by the MARCS model is significantly smaller than observed. This is because the 
MARCS model with the parameters of Aldebaran has a geometrical extension of only ~2% (with respect to the stellar radius). The 
observed spectrum and interferometric data in the CO lines can be simultaneously reproduced by placing an additional CO layer 
above the MARCS photosphere. This CO layer is extended to 2.5 ± 0.3 if* with CO column densities of 5 x 10 19 -2 x 10 20 cnr 2 and 
a temperature of 1500 ± 200 K. 

Conclusions. The high spectral resolution of AMBER has enabled us to spatially resolve the inhomogeneous, extended outer atmo- 
sphere (MOLsphere) in the individual CO lines for the first time in a K giant. Our modeling of the MOLsphere of Aldebaran suggests 
a rather small gradient in the temperature distribution above the photosphere up to 2-3 R+. 



Key words, infrared: stars 
Aldebaran 



techniques: interferometric - stars: mass-loss - stars: late-type - stars: atmospheres - stars: individual: 



> 

m 

rn 
O 
m 



x 



1. Introduction 

Red giant stars, which represent the late evolutionary stages of 
intermediate- to low-mass stars from the red giant branch (RGB) 
to the asymptotic giant branch (AGB), experience mass loss 
ranging from 10~" Mq yr~' up to 10~ 4 Mq yr . However, the 
mass loss in red giants is not yet understood well. For Mira-type 
AGB stars with large variability amplitudes (AV « 9), the com- 
bination of the levitation of the atmosphere by the stellar pulsa- 
tion and the radiation pressure on dust grains is often believed to 
drive the mass loss (see, e.g., Hofner l201 II for discussion of the 
success and the current problem with this mechanism). However, 
it is by no means clear whether this mechanism can operate in 
red giants in general, that is, normal (i.e., non-Mira-type) K and 
M giants, whose variability amplitudes (AV ~ 1-2) are signif- 
icantly smaller than Miras. Furthermore, K giants and early M 
giants are experiencing mass loss, despite the absence of dust. 

Observations of the outer atmosphere are important for un- 
derstanding the mass-loss mechanism. Infrared spectroscopic 
studies of molecular lines reveal the presence of an extended 
molecular outer atmosphere, the so-called MOLsphere, in K 
and M giants, as well as in red supergiants (Tsuji 1988. 2000a 
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2000b, 2001; Tsuji et al. [l997l Spectro-interferometry, which 
combines high spectral resolution and high spatial resolution, is 
powerful for studying the physical properties of the MOLsphere. 

We spatially resolved the normal M7 giant BK Vir in the 
CO first overtone lines near 2.3 yum with a spectral resolution of 
12000 using VLTI/AMBER (Ohnaka et al. [20T2b . The observed 
CO line spectrum can be reproduced well by the MARCS photo- 
spheric models (Gustafsson et al. 2008), which are state-of-the- 
art models for the photosphere of cool evolved stars. However, 
the MARCS models clearly fail to explain the observed angular 
size of the star in the individual CO lines. The observed angular 
size in the CO lines is much greater than the models predict, sug- 
gesting the presence of more extended components. We found 
out that the observed spectrum and angular size can be explained 
by a model in which two extra CO layers are added above 
the MARCS model photospheres, at 1.2-1.25 tf* and 2.5-3 fl* 
with CO column densities of 10 22 ctrT 2 and 10 19 — 10 20 cm -2 . 
Surprisingly, the temperatures of the CO layers (1900-2100 K 
and 1500-2100 K) are higher than or equal to the temperature of 
the uppermost layer of the photosphere (see Fig. 9 of Ohnaka et 
al. 120121 1. While many K and M giants possess chromospheres, 
the derived temperatures of the MOLsphere are still lower than 
in the chromosphere (-6000-10000 K). Therefore, although the 
MOLsphere shows a temperature inversion, it is still different 
from the chromosphere. 
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Table 1. Summary of the AMBER observation of Aldebaran with the B2-C1-D0 AT configuration. The seeing and coherence time 
(tq) are in the visible. 
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121 


2500 
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Fig. 1. AMBER observation of Aldebaran in 
the CO first overtone lines, a: Visibility ob- 
served on the 10.4 m baseline (red solid line), 
b: Uniform-disk (UD) diameter (red solid line), 
c: Differential phase (red solid line). In each 
panel, the scaled observed spectrum is shown 
by the black solid line. The positions of the CO 
lines are also marked with the ticks. 



A better understanding of the physical properties of the 
MOLsphere is important for clarifying the mass-loss mecha- 
nism, because the physical process responsible for forming the 
MOLsphere may be related to the driving mechanism of the 
mass loss. Given that IR spectro-interferometric observations of 
normal K and M giants are still scarce, AMBER observations of 
more normal red giants are indispensable for shedding new light 
on the origin of the MOLsphere. 

As a first step, we present a high-spatial and high-spectral 
resolution observation of the CO first overtone lines in the well- 
studied K5 giant Aldebaran (a Tau) with VLTI/ AMBER. At 
a distance of 20.4 pc (parallax = 48.94 mas, van Leeuwen 
120071 1. Aldebaran is one of the nearby red giants. Its angular 
diameter has been measured by the lunar occultation and long- 
baseline optical/IR interferometry. As summarized in Richichi 
& Roccatagliata (2005), most of the angular diameter measure- 
ments between the V and L band indicate 18-22 mas but without 
a noticeable wavelength dependence. As discussed in their pa- 
per, the scatter in the measured angular diameters mostly reflects 
the discrepancies among the lunar occultation measurements, 
and the scatter among the diameters measured by long-baseline 
interferometry is smaller. In the present work, we adopt the 
uniform-disk diameter of 19.96+0.03 mas derived by Richichi & 
Roccatagliata (2005), which is the average of their lunar occul- 
tation measurements at 2.22 and 3.55 //m and the VLTI/VINCI 
measurements in the K band. They argue that their measure- 
ments would represent the best available value for the angular 
diameter of Aldebaran. We adopt an effective temperature (T e g) 



of 3874 + 100 K and a stellar mass (M*) = 1.5 + 0.3 M Q from 
Tsuji (2008 ) for the present work. Combining this stellar mass 
and the radius of 44 Rq measured by Richichi & Roccatagliata 
(2005) results in a surface gravity of log g =1.5 (in units of 
cm s~ 2 ). The metallicity of Aldebaran is solar or marginally sub- 
solar with [F e/H] = -0. 14 + 0.30 (Kovacs fl983l) and -0. 15 + 0.2 
(Decin et al.|2003 ). Based on high-resolution spectra in the H, 
K, and L bands, Tsuji (2008 ) derived the carbon, nitrogen, and 
oxygen abundances to be logAc = 8.38 + 0.04, logAN = 8.05, 
and log Aq = 8.79 + 0.04 (abundances are expressed on the loga- 
rithmic scale with log Ah = 12). These CNO abundances suggest 
the mixing of the CN-cycled material in the first dredge-up, al- 
though the derived 12 C/ 13 C ratio of 10 cannot be fully understood 
by the standard evolutionary theory. The micro-turbulent veloc- 
ity derived in previous studies is approximately 2 km s -1 (Smith 
& Lambert [191351 Decin et al. l200ll Tsuii l2008l . The mass-loss 
rate and the terminal wind velocity of Aldebaran are estimated 
tobe (l-1.6)xl0' 11 M Q yr" 1 (Robinson et al. [19981 Wood et al. 
I2007T) and 30 km s _1 (Robinson et al. [1998b . respectively. 

Tsuji (1200 11 1 discovered the H2O absorption feature at 
~6.6 fim in Aldebaran, which originates in the MOLsphere. 
Based on a detailed analysis of a high-resolution spectrum of 
Aldebaran, Tsuji (2008 ) revealed absorption excess in the strong 
CO first overtone lines, which is also interpreted as the contribu- 
tion of the MOLsphere. On the other hand, the emission lines in 
the UV (e.g., Mgn, Fen, and Oi lines) indicate the presence of 
the chromosphere. This inhomogeneous, multi-component na- 
ture of the outer atmosphere is consistent with what is suggested 
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from the observations of the CO fundamental lines in the IR 
(Wiedemann et al. 1994) and the CO fluorescent lines in the UV 
(McMurry & Jordan 2000 ). Note, however, that the nature of the 
CO gas probed by the 4.6 //m fundamental lines (v = 1-0) 
might be different from that probed by the first overtone lines 
(u — 2 - 0) near 2.3 fim in the present work. The spectral energy 
distribution (SED) of Aldebaran from the optical to ~1 mm, as 
well as the mid-IR spectrum shows no signature of dust (Dehaes 
et al. 1201 It Monnier et al. 11 9981 ), and therefore, this star is ideal 
for studying the mass-loss mechanism without dust. 

2. Observation 

AMBER is a spectro-interferometric instrument at VLTI and can 
provide high spatial resolutions and high spectral resolutions of 
up to 12000 (Petrov et al. 120071 ) by combining three 8.2 m Unit 
Telescopes (UTs) or 1.8 m Auxiliary Telescopes (ATs). AMBER 
measures the visibility, which is the amplitude of the Fourier 
transform of the object's intensity distribution on the sky, as 
well as the closure phase and differential phase. The closure 
phase contains information about the asymmetry of the object, 
with non-zero and non-7r closure phases being the signature of 
asymmetry. The differential phase provides information about 
the photocenter shift of the object in spectral features with re- 
spect to the continuum. 

We observed Aldebaran on 2012 December 12 with AMBER 
using the B2-C1-D0 (11— 23-34m) linear array (Program ID: 
090.D-0459A, PL: K. Ohnaka), as summarized in Table [Q We 
used the high-spectral resolution mode in the K band (HR_K) 
with a spectral resolution of 12000 between 2.28 and 2.31 fim 
to observe the CO first overtone lines near the 2-0 band head 
at 2.294 fim. The VLTI fringe tracker FINITO was not used, be- 
cause while the visibility in the H band, where FINITO operates, 
is higher than the limit of 0.15 on the shortest baseline, the H- 
band visibility on the second shortest baseline is below the limit. 
However, the high brightness of Aldebaran enabled us to detect 
fringes without FINITO with a Detector Integration Time (DIT) 
of 121 ms. We observed Procyon (F5IV-V) not only as the inter- 
ferometric calibrator but also as the spectroscopic calibrator. 

We reduced the data using amdlib ver.3.0.53 as described in 
Ohnaka et al. (2009]). We took the be st 80% of the frames in 
terms of the fringe S/N (Tatulli et al. 120071 ) . For the calibra- 
tor Procyon, we adopted the uniform-disk diameter of 5.38 + 
0.05 mas from the CHARM2 catalog (Richichi et al. 120051 
There are two issues that had to be taken into account in the 
reduction of the present data. Firstly, the telescope at the DO sta- 
tion was shadowed by UT1 during the observation of Aldebaran 
(but not during the observation of Procyon). Therefore, the flux 
of Aldebaran at this telescope was lower than at the other two 
telescopes by a factor of -10, which makes the derivation of the 
interferometric observables on the baselines including the DO 
telescope less reliable. Therefore, in the present work, we only 
discuss the visibility and differential phase obtained on the B2- 
Cl baseline, which provides a projected baseline length (B p ) of 
10.4 m and a position angle (PA) of 37°. Secondly, only one 
data set of the calibrator was obtained, which makes it difficult 
to carry out a reliable calibration of the absolute visibility level. 
Because the time variation expected in the angular diameter of 
Aldebaran is as small as +0.17 mas (Richichi & Roccatagliata 
2005), we scaled the observed calibrated visibility in the con- 
tinuum to those expected from a uniform disk with the adopted 
angular diameter of 19.96 mas. 

1 Available at http://www.jmmc.fr/data_processing_amber.htm 



The wavelength calibration was carried out using the tel- 
luric lines identified in the observed spectrum of Procyon, as 
described in Ohnaka et al. (120121 ). The uncertainty in the wave- 
length calibration is 3.2 X 10~ 5 fim (4.2 km s~'). The wavelength 
scale was converted to the laboratory frame using a heliocentric 
velocity of 54.2 km s of Aldebaran (Gontcharov 2006). The 
spectrum of Procyon in the observed spectral window is feature- 
less with AMBER's spectral resolution except for the Mg line at 
2.2814 //m. Therefore, the calibrated spectrum of Aldebaran was 
obtained by dividing the observed spectrum of Aldebaran with 
that of Procyon. 



3. Results 

Figure QJ shows the observed visibility in the CO lines on the 
10.4 m baseline with the observed (spatially unresolved) spec- 
trum. The signature of the CO lines is clear in the observed 
visibility. The uniform-disk diameter derived from the visibil- 
ity is shown in Fig. QJ). The uniform-disk diameter in the CO 
lines is 20-35% larger than in the continuum. Figure [TJ: shows 
the observed differential phase. We detect non-zero differential 
phases in the CO lines (i.e., wavelength-dependent photocenter 
shifts), which indicate the presence of asymmetric or inhomo- 
geneous structures in the CO-line-forming region on a spatial 
scale of ~45 mas (corresponding to the spatial resolution with 
the 10.4 m baseline). This is more than twice as large as the 
angular diameter of the star itself and consistent with the size of 
the MOLsphere as discussed in Sect. [4] Richichi & Roccatagliata 
(2005) found possible surface structures — though not conclusive 
as the authors mention — in the 1-D image of Aldebaran recon- 
structed from the lunar occultation data taken at 2.22 //m with a 
spectral resolution of ~70, which sample the spectral region free 
from strong lines. Our AMBER observation is the first study to 
spatially resolve the atmosphere of a K giant in individual CO 
lines and detect asymmetry in the CO-line-forming region. 

The observed increase in the uniform-disk diameter in the 
CO lines is in marked contrast to the negative detection of 
a change in the angular size in the TiO band at 7120 A 
(Quirrenbach et al. 1993 ). However, this might be due to the 
lower spectral resolution used in their observations. For exam- 
ple, if we bin our AMBER data to a spectral resolution of 1500 
as described in Ohnaka et al. (2009 ), the increase in the uniform- 
disk diameter in the CO band head is 4%. The narrow-band filter 
for the TiO band used by Quirrenbach et al. (11993) has a spec- 
tral resolution of ~60. The width of the filter (120 A) is smaller 
than the width of the TiO band (-400 A) starting at -7050 A. 
However, the observed spectrum of Aldebaran in Kieling (1987 
Fig. 4c, HR 1457 = Aldebaran) shows that this TiO band is cen- 
tered at 7200 A, and the central wavelength of the filter (7 120 A) 
does not exactly match the center of the TiO band. The filter 
samples the spectral region between the beginning and the cen- 
ter of the TiO band, not just the deepest part of the TiO band. 
Therefore, the low spectral resolution, as well as the shift be- 
tween the filter's central wavelength and the center of the TiO 
band, might have led to the negative detection of the extended 
atmosphere of Aldebaran. It might not have been a problem for 
detecting a change in the angular size for cooler stars with more 
pronounced extended atmospheres. In any case, measurements 
of the angular diameter of K giants in the TiO bands with higher 
spectral resolutions are necessary to compare with the angular 
diameter in the CO first overtone lines. This can be feasible, for 
example, with the visible interferometric instrument VEGA at 
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CHARA with spectral resolutions of up to 30000 (Mourard et 
al. 12009b . 

The high spatial and high spectral resolution of AMBER has 
enabled Ohnaka et al. (20091 |20TT1 and l20T3b to spatially re- 
solve the gas motions in the photosphere and MOLsphere in the 
red supergiants Betelgeuse and Antares. However it is difficult 
to carry out a similar analysis for the current data of Aldebaran, 
because the amplitude of the velocity field in the photosphere 
(and possibly also in the MOLsphere) is much smaller than in 
the red supergiants. For example, Gray (2009 1 measured a time 
variation in the radial velocity of +0.3 km s _I in Aldebaran, 
which is much smaller than the amplitude of ~7 km s _1 found 
in Betelgeuse (Gray 2008). The micro- and macro-turbulent ve- 
locities of ~2 km s _1 and 3.6 + 0.3 km s _1 derived for Aldebaran 
(references given in Sect. [T]and Tsuji 1986) are also noticeably 
smaller than the micro- and macro-turbulent velocities of ~5 
and 10-20 km s" 1 in Betel geuse (Ohnaka et al. 120091 and ref- 
erences therein). The amplitude of the velocity field suggested 
from these observations is too small compared to AMBER's 
spectral resolution and the accuracy of the wavelength calibra- 
tion. 



4. Modeling the AMBER data 

We first compare the observed spectrum and visibility with 
the MARCS photospheric model (Gustafsson et al. 12008b . The 
MARCS models represent plane-parallel or spherical photo- 
spheres in hydrostatic and radiative equilibrium, with the effects 
of molecular and atomic line opacities included using the opacity 
sampling method. Each MARCS model is specified by effective 
temperature (T e ff), surface gravity (log g), micro-turbulent veloc- 
ity (fmicroX chemical composition, and (for the spherical case) 
stellar mass (M*). We selected spherical MARCS models with 
T eff = 3900 K, log g = 1.5, M* = 1 and 2 M Q (no MARCS model 
with 1.5 M* is available), i> m j cro = 2 km s _1 , and the "moderately 
CN-cycled" composition with [Fe/H] = 0.0, because these mod- 
els have the parameters closest to those of Aldebaran summa- 
rized in Sect. Q] 

Using the temperature and pressure distributions from the 
MARCS models, we computed the monochromatic intensity 
profile and visibility, as well as the spectrum using the CO line 
list of Goorvitch (1994), as described in Appendix lAl The an- 
gular scale of the models was determined so that the uniform- 
disk diameter derived from the model visibility in the continuum 
agrees with the adopted value of 19.96 mas. 

Comparison of the MARCS model visibility and spectrum 
with the AMBER observation is shown in Fig. [2] The observed 
CO line spectrum (Fig. |2^) is reproduced well by the MARCS 
model with T eS = 3900 K, logg = 1.5, and M* = 1 M Q . 
However, as Figs. \2)p and [2J; show, the model predicts only 
marginal decrease in the visibility (thus only marginal increase 
in the uniform-disk diameter) in the CO lines and fails to ex- 
plain the observed data. This means that the star is much more 
extended than the MARCS model predicts. The geometrical ex- 
tension of the MARCS model, which is defined as the geometri- 
cal distance between the layers with tr oss = 1 and tr oss = 10~ 5 
(tross is the optical depth based on the Rosseland mean opacity), 
is 2% of the stellar radius. This geometrically thin photosphere 
is seen in Fig. [2jl, where the intensity profiles in the continuum 
and in the CO band head are shown. This means that the visibil- 
ities predicted in the continuum and in the CO band head differ 
only slightly (Fig. |2^). This is why the MARCS model cannot 
explain the increase in the uniform-disk diameter of 20-35% in 
the CO lines. We also examined whether MARCS models with 
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Fig. 2. Comparison of the MARCS model with the AMBER 
data of Aldebaran. The parameters of the MARCS model are 
given in Sect. [4] a: Model spectrum (black solid line) computed 
from the pressure and temperature distributions of the MARCS 
model as described in Appendix [A] (not the pre-computed spec- 
trum available on the MARCS website) with the observed spec- 
trum (red solid line). The green and blue dots show the wave- 
lengths of the intensity profiles and visibilities shown in the pan- 
els d and e. The positions of the CO lines are marked by the 
ticks, b: Model visibility predicted at the 10.4 m baseline (black 
solid line) with the observed data (red solid line), c: Uniform- 
disk (UD) diameter derived from the MARCS model visibility 
shown in the panel b (black solid line) with the one observed 
(red solid line), d: Model intensity profiles predicted in the con- 
tinuum (green solid line) and in the CO band head (blue solid 
line), at the wavelengths marked in the panel a. e: Model visi- 
bilities as a function of spatial frequency predicted in the contin- 
uum (green solid line) and in the CO band head (blue solid line). 
The visibilities observed at these wavelengths are plotted by the 
dot (continuum) and the filled diamond (CO band head). 



slightly different parameters can explain the observed visibility. 
However, the MARCS models with T eS = 3800^1000 K, log^ = 
1.0 and 1.5, and M* = 1 and 2 Mq are characterized by the geo- 
metrical extension of 4% at most, and therefore, cannot explain 
the observed significant increase in the angular size in the CO 
lines. This finding that the MARCS model can reproduce the CO 
first overtone line spectrum well but fails to explain the angular 
size of the star in the CO lines is the same as what we found in 
the M7 giant BK Vir in Ohnaka et al. ( f2012b . Our AMBER ob- 
servation with high spectral resolution has clearly revealed the 
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presence of the MOLsphere, which cannot be accounted for by 
the photospheric models. 

We derived the physical properties of the MOLsphere us- 
ing the model described in Ohnaka et al. (120121 1. In this model, 
one or two CO layers are added above the MARCS photospheric 
model. Given that the visibility was obtained only on one base- 
line in the present work, we used only one layer to keep the 
number of free parameters as small as possible. The geometrical 
thickness of the additional CO layer is fixed to 0.1 and the 
micro-turbulent velocity in the CO layer is set to 2 km s _1 , which 
is the same as in the MARCS photospheric model. The (inner) 
radius, CO column density in the radial direction, and tempera- 
ture of the CO layer were treated as free parameters. The com- 
putation of the intensity profile, visibility, and spectrum from the 
MARCS +MOLsphere model is described in Appendix lAl 

Figure [3] shows a comparison of the best-fit 
MARCS +MOLsphere model with the observed spectrum 
and visibility. The MOLsphere of this model is characterized 
by a radius of 2.5 with a temperature of 1500 K and a CO 
column density of 1 x 10 20 crrT 2 . The observed spectrum and 
visibility, as well as the uniform-disk diameter are reasonably 
reproduced. The intensity profile shown in Fig. [3jl reveals the 
extended emission from the MOLsphere. While the intensity 
of this extended emission is low compared to the intensity on 
the stellar disk, the flux contribution is not negligible because 
it is very extended. As a result, the MOLsphere makes the 
star appear much more extended than the photosphere alone 
and can explain the observed angular size in the CO lines. 
We computed MOLsphere models with radii between 1.2 and 
2.8 R+, temperatures between 1000 and 2000 K, and CO column 
densities between 1 x 10 19 and 1 x 10 21 cm~ 2 , and estimate the 
uncertainties in the radius, temperature, and CO column density 
to be +0.3 R+, +200 K, and a factor of 2, respectively. 

It may appear contradictory that the model spectrum is 
hardly affected by the MOLsphere, while the angular size is 
largely affected. This can be explained as follows. Figure [4] 
shows spatially resolved spectra predicted at two positions (at 
the center of the stellar disk and at the edge of the MOLsphere) 
by the MARCS+MOLsphere model, together with the spatially 
unresolved model spectrum. On the one hand, the additional CO 
layer introduces stronger absorption (i.e., CO absorption lines 
are deeper) for lines of sight within the stellar disk. This can be 
seen in the spatially resolved spectrum at the stellar disk center 
(red line) when compared to the spatially unresolved spectrum 
from the MARCS-only model (black line). On the other hand, 
the spatially resolved spectrum at the edge of the MOLsphere 
(blue line) shows that the CO lines appear in emission off the 
limb of the star (note that the off-limb spectrum is scaled up 
by a factor of 10 for visual clarity). This extended emission 
from the MOLsphere fills in the additional absorption due to the 
MOLsphere itself, which results in little change in the spatially 
unresolved spectrum. The spectra predicted by the MARCS- 
only and MARCS+MOLsphere models shown in Figs. [2^ and 
[3^ show very little difference. However, the difference in the 
visibility predicted by these models is obvious. This illustrates 
that the MOLsphere can become invisible for a spectrometer 
but not for a spectro-interferometer (however, signatures of the 
MOLsphere can be detected by spectroscopy for some lines, as 
demonstrated by Tsuii l200"Tl [2008). It is also feasible to observe 
the CO emission line spectrum expected from the MOLsphere, if 
we carry out aperture-synthesis imaging with high-spectral res- 
olution with AMBER, as demonstrated for the red supergiants 
Betelgeuse and Antares (Ohnaka l2013l Ohnaka et al. 1201 31 >. 
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Fig. 3. Comparison of the MARCS+MOLsphere model with 
the AMBER data of Aldebaran shown in the same manner as 
in Fig. [2] 



We assumed local thermodynamical equilibrium (LTE) for 
modeling the MOLsphere. We examined whether this assump- 
tion is valid for the derived density and temperature of the 
MOLsphere. As described in Ohnaka (2004), we compared the 
collisional and radiative de-excitation rates. The collisional de- 
excitation rate C u \ is estimated as Ncr u \v K i, where N, cr ul , and u re i 
denote the number density of the primary collision partner, col- 
lisional cross section, and relative velocity between the collision 
partner and CO molecules. We approximate <x ul with the geomet- 
rical cross section of 10~ 15 cm 2 and adopt u,ei = 5 km s _1 . The 
number density of the primary collision parter (H2 or H) is esti- 
mated as follows. We estimate the CO number density by divid- 
ing the derived CO column density by the geometrical thickness 
of the MOLsphere. However, we arbitrarily set the geometrical 
thickness to 0. 1 R* because this parameter cannot be constrained 
from the current data. Therefore, we estimated a lower limit on 
the CO number density by dividing the derived CO column den- 
sity (10 20 ctrT 2 ) by the geometrical distance between the stel- 
lar surface and the radius of the MOLsphere. This results in a 
CO number density of 2.2 x 10 7 cm 4 . In chemical equilibrium 
at 1500 K, this CO number density suggests H2 and H number 
densities of 5.6 x 10 10 and 9.5 x 10 9 cm , respectively. Using 
the H2 number density for N, we obtain C u \ = 2.8 x 10 s . This 
collisional de-excitation rate is much higher than the radiative 
de-excitation rates A u \ - 0.6-1 s _1 . Therefore, it is unlikely that 
our modeling is significantly affected by NLTE effects. 
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Fig. 4. a: The image predicted at the CO band head by the 
MARCS +MOLsphere model shown in Fig. [3] b: Spatially re- 
solved spectra at the center of the stellar disk (position A in 
the panel a) and at the edge of the MOLsphere (position B) are 
shown by the red and blue solid line, respectively. The spatially 
unresolved spectrum from the MARCS -only model is plotted by 
the black solid line. The spectrum at the edge of the MOLsphere 
is scaled up by a factor of 10. 



5. Discussion 

Our modeling of the spatially and spectrally resolved CO line 
data suggests a high CO column density in the MOLsphere of 
Aldebaran. Tsuji ((2001) estimated an H2O column density of 
2 x 10 17 cm" 2 from the weak absorption feature near 6.6 //m. The 
CO number density estimated above suggests an H2O column 
density of 1.4 x 10 20 cm" 2 , if we assume chemical equilibrium 
at 1500 K. The H2O column density implied from our modeling 
of the CO lines is much higher than estimated from the previ- 
ous spectroscopic study. A possible reason for this disagreement 
may be that it is not straightforward to disentangle the effects of 
the geometrical extension and column density from the observed 
spectra alone. 

Based on the spectral analysis of the CO fluorescent lines in 
the UV, McMurry & Jordan (2000 ) deduced the inhomogeneous 
structure of the atmosphere of Aldebaran: the chromosphere as 
hot as 10 4 K extending into the transition region with 10 s K and 
the cool CO gas at ~2000 K. This temperature of the cool CO 
gas is comparable to the temperature of the MOLsphere derived 
above. However, the CO column density derived by McMurry & 



10" 



5000 



2000 



1000 



4000 



3000 



-2000 



1000 



I MARCS photosphere 



MOLsphere 



radius(i?J 



MARCS photosphere 



MOLsphere 



1 .5 



2 

radius(Rj 



2.5 



Fig. 5. Temperature distribution as a function of radius predicted 
by the MARCS model for Aldebaran (black solid line) and the 
temperature of the MOLsphere derived from our modeling (red 
hatched rectangle). The panel a shows the entire temperature 
range covered by the MARCS photosphere (note the logarith- 
mic scale both for the temperature and radius), while the panel 
b shows a lower temperature region including the upper photo- 
sphere and MOLsphere on the linear scale. 



Jordan (2000) is ~ 1 x 10 17 cm" 2 , which is much lower than de- 
rived for the MOLsphere. We computed MARCS +MOLsphere 
models with the CO column density fixed to 1 x 10 17 cm" 2 , but 
the angular size in the CO lines predicted by these models is too 
small compared to the observed data, simply because the den- 
sity is too low. Therefore, this discrepancy in the CO column 
density cannot be reconciled by adjusting the parameters in our 
modeling of the MOLsphere. Instead it may be related to the 
location where the CO fluorescent lines originate. The CO fluo- 
rescent lines in the UV are pumped by the 1 lines that form in 
the chromosphere at temperatures of 6500-8000 K (McMurry 
1999; McMurry & Jordan l2000l . The spectral analysis of the 
latter authors implies that "the 1 radiation and CO molecules 
must be in close proximity". Therefore, the CO column density 
derived from the fluorescent lines may refer only to the cool CO 
gas near the hot chromospheric gas, while the CO column den- 
sity estimated from the CO first overtone lines refers to the entire 
cool CO gas. 

Figure [5] shows the temperature distribution of the MARCS 
model for Aldebaran and the temperature and location of the 
MOLsphere derived from our modeling. Compared to the steep 
temperature decrease in the photosphere, the temperature gra- 
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dient in the outer atmosphere is much smaller. This is similar 
to what was found in the M7 giant BK Vir, in which the tem- 
perature in the MOLsphere is equal to or even higher than in the 
uppermost layer of the MARCS model (Ohnaka et al. 120121 . The 
CO column density in the MOLsphere of Aldebaran is compara- 
ble to the value found in the red supergiants Betelgeuse (Ohnaka 
et al. 120091 [20TTb and Antare s (Ohnaka et al. |2DT3T >. as well as 
BK Vir (Ohnaka et al. 2012), despite the higher effective tem- 
perature and higher surface gravity of Aldebaran. This implies 
that the physical process responsible for the formation of the 
MOLsphere may be the same for red supergiants and (non-Mira) 
red giants regardless the spectral type and luminosity class. 

To shed new light on the origin of the MOLsphere, the mag- 
netohydrodynamical (MHD) simulation of Suzuki (2007) is in- 
teresting for the present study. In particular, his model IV has 
stellar parameters similar to Aldebaran: r e s = 3900 K, M* = 
1 Mq, R+ = 3 1 Rq, and log g — 1.4. He simulates Alfven waves, 
which are excited by the surface convection and travel outward 
along an open 1-D magnetic flux tube from the photosphere 
to ~25 The simulation also solves the thermal structure of 
the wind, instead of assuming a given temperature profile. The 
magnetic field used in the simulation corresponds to a surface- 
averaged magnetic field strength of 1 G. The stellar wind of 
the model IV shows hot gas bubbles with a temperature higher 
than ~ 10 5 K embedded in much cooler gas. A snapshot of the 
temperature structure (see Fig. 9 in his paper) shows that the 
temperature in the cool gas can be as low as 1000-2000 K, 
which agrees with the observationally derived temperature of 
the MOLsphere. The model also show hot gas bubbles with tem- 
peratures of 1 4 — 1 5 K, which may correspond to the chromo- 
sphere or the "buried" corona as Ayres et al. (2003 1 propose. 
The lower limit on the CO number density of 2.2 x 10 7 cm -3 
derived above suggests a gas density of 1.9 x 10~ 13 g cm -3 . The 
time-averaged gas density predicted at 2 for the model IV is 
~10~ 14 g cm 4 . Given the large temporal fluctuation in the den- 
sity structure predicted by the model (by a factor of 10-100) and 
the uncertainty in the observationally derived CO column den- 
sity, this value roughly agrees with the observationally estimated 
density. 

However, the time-averaged mass-loss rate predicted by this 
model is ~ 10~ 9 Mq yr~', which is significantly higher than 
the observationally estimated value of (1-I.6)xl0~ n Mq yr _1 
(Robinson et al. 119981 Wood et al. l2007l l. although the model is 
not constructed particularly for Aldebaran. Also, as pointed out 
by Airapetian et al. (2010), the fully ionized plasma is assumed, 
while the winds from red giants are weakly ionized. While the 
MHD simulation of the stellar wind of Aldebaran by Airapetian 
et al. (120101 ) assumes weakly ionized plasma, the temperature 
structure is not self-consistently computed in their simulation, 
and therefore, their model cannot yet account for the MOLsphere 
component. Another uncertainty in the MHD models is that the 
strength of the magnetic fields in Aldebaran is not observation- 
ally determined. On the one hand, the detection of O vi emission 
lines suggests that the outer atmosphere of Aldebaran is influ- 
enced by magnetic fields (Harper et al. 1201 II . On the other hand, 
the spectropolarimetric observations of Auriere et al. (2010jl de- 
tected no signature of magnetic fields in Aldebaran, implying 
that the surface-averaged field strength is much smaller than the 
1 G in the model. An extension of the MHD simulations of 
Suzuki (120071 1 with weaker magnetic fields and weakly ionized 
plasma, as well as the inclusion of the MOLsphere component in 
the model of Airapetian et al. d2010l l would be crucial for com- 
parison with the observed data. 



McMurry & Jordan (120001 propose an alternative scenario, 
in which the hot chromospheric gas is associated to acoustic 
shocks. The one-dimensional simulation of acoustic waves for a 
K5 giant (corresponding to Aldebaran) by Buchholz et al. (1998 ) 
shows that while the temperature behind shocks can become 
higher than -5000 K, the gas between shocks can be as cool 
as 2000 K. However, the cooling due to molecules (most impor- 
tantly CO) is not included in their simulation. A more realistic 
simulation of acoustic waves including the molecular cooling is 
necessary to understand the formation mechanism of the chro- 
mosphere and MOLsphere. 

While our MARCS +MOLsphere model simultaneously ex- 
plains the spectrum and angular size in the CO lines of 
Aldebaran, there is still a problem with the MOLsphere as men- 
tioned by Ohnaka et al. (|2012] | for the M7 giant BK Vir. A 
large number of other molecular species is expected from the 
MOLsphere, most notably TiO. The contribution of TiO in the 
MOLsphere can make the TiO absorption bands too strong com- 
pared to the observed spectra. Ohnaka et al. (2012) suggest that 
the TiO bands may form by scattering in the MOLsphere, and 
the scattered emission may fill in the absorption caused by the 
MOLsphere. However, it is necessary to construct a MOLsphere 
model with the scattering taken into account and compare it with 
the observed spectra and the available and future angular diam- 
eter measurements in the TiO bands. 



6. Conclusions 

We have spatially resolved the K5 giant Aldebaran in the in- 
dividual CO first overtone lines near 2.3 pm, taking advantage 
of high spatial and high spectral resolution of VLTI/AMBER. 
The observed uniform-disk diameter in the CO lines is 20-35% 
larger than in the continuum. This increase in the angular size 
in the CO features would be difficult to detect with a spectral 
resolution lower than ~1500 (AMBER's medium spectral reso- 
lution mode). Our observation of Aldebaran illustrates the power 
of the high-spectral resolution capability of AMBER. While the 
MARCS photospheric model can reproduce the observed CO 
line spectrum very well, it fails to explain the observed angular 
diameter in the CO lines: the model predicts the angular diam- 
eter to be too small compared to the observation. This reveals 
the presence of the MOLsphere in a K giant for the first time in 
the individual CO lines. Our modeling of the AMBER data sug- 
gests that the MOLsphere extends to ~2.5 R* with a temperature 
of -1500 K and a CO column density of ~1 x 10 20 cm" 2 . The 
temperature gradient in the MOLsphere is much smaller than 
in the photosphere. The derived temperature and density of the 
MOLsphere roughly agree with the MHD simulation of Suzuki 
(2007), although a modeling with more appropriate assumptions 
for Aldebaran is necessary for understanding the origin of the 
MOLsphere. 

Our work demonstrates that the high spatial and high spec- 
tral resolution of VLTI/AMBER enables us to extract unique in- 
formation on the outer atmosphere even from data on only one 
baseline. Now it is crucial to carry out a systematic survey of 
red (super)giants with a range of effective temperature and lumi- 
nosity with this technique for obtaining a more comprehensive 
picture of the dependence (or absence of it) of the properties of 
the MOLsphere on stellar parameters. This is possible with the 
current AMBER's sensitivity. 
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where S a and ta denote the source function and the optical depth 
at A along the ray with the impact parameter p, and the inte- 
gration is carried out along this ray. In local thermodynamical 
equilibrium (LTE), the source function is given by 



K.l 



o-.i 



KA +<Ta 



-./.I 



where ka and cr A represent the absorption and scattering coeffi- 
cients, respectively. B \ and J a represent the Planck function and 
the mean intensity, respectively. From the pressure and temper- 
ature distributions of a given MARCS model, we compute the 
number of each molecular and atomic species in chemical equi- 
librium. With the Voigt function adopted for the line profile, the 
absorption coefficient at each wavelength is calculated from the 
CO line list of Goorvitch (1994) and the continuous opacities 
due to the bound-free and free-free transitions of H and H, the 
free-free transitions H, 1 and He -1 , and the bound-free transi- 
tions of Si, Mg, and Ca. We used the data of John ( 119881 ) for the 
H" 1 opacity, while we adopted the cross sections given in Tsuji 
(|1971| and references therein) for the other continuous opacities. 

The intensity profile from a MARCS +MOLsphere model 
4>iarcs+mol0?, X) is computed using the MARCS-only intensity 
profile /marcsCp, X) as follows: 

^marcs+mol(/?, J) — ^marcs(p^) e ^ + ^(TmolXI - e Tj ), 

where r^ OL denotes the optical depth in the MOLsphere along 
the ray with the impact parameter p, and Tmol is the temperature 
of the MOLsphere. 

Once the intensity profile at each wavelength is obtained, the 
spectrum from the MARCS +MOLsphere model (Fa) is calcu- 
lated by 



■x 



Imarcs+mol(p, T) M dp., 



where p is defined as y/l - (p/Rmol) 2 with Rmol being the 
outer radius of the MOLsphere. The visibility is obtained by tak- 
ing the Hankel transform (2-D Fourier transform for axisymmet- 
ric objects) of the intensity profile. 

Finally, the monochromatic intensity profile and visibility, as 
well as the spectrum are spectrally convolved with the AMBER's 
spectral resolution. 



Appendix A: Computation of the intensity profile, 
spectrum, and visibility from the 
MARCS+MOLsphere model 

We first compute the monochromatic intensity profile only from 
a given MARCS model with a wavelength interval of 2.7 x 
10~ 6 pm. This wavelength interval was chosen to sample at 
least 4-5 points across the line profile whose width is deter- 
mined by the thermal and micro-turbulent velocity. As depicted 
in Fig. IA.1I the monochromatic intensity 7marcs(p, A) at the 
wavelength A and at a given impact parameter (p) is obtained 
by 



4iarcs(/?,T)= S a (ta) e TA dr A , 
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Fig.A.l. Schematic view of the computation of the intensity 
profile from the MARCS +MOLsphere model. Each MARCS 
model consists of 56 layers, but only six layers are drawn for 
the sake of clarity. 



